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Ground-state hyperfine splittings in hydrogen and muonium are very well measured. Their difference, after 
correcting for magnetic moment and reduced mass effects, is due solely to proton structure — the large QED con- 
tributions for a pointlike nucleus essentially cancel. The rescaled hyperfine difference depends on the Zemach 
radius, a fundamental measure of the proton, computed as an integral over a product of electric and magnetic 
proton form factors. The determination of the Zemach radius, (1.043 ± 0.016) fm, from atomic physics tightly 
constrains fits to accelerator measurements of proton form factors. Conversely, we can use muonium data to 
extract an "experimental" value for QED corrections to hydrogenic hyperfine data; we find that measurement 
and theory are consistent. 



Introduction. Quantum Electrodynamics, QED, stands out 
as the most precisely tested component of the Standard Model. 
QED predictions for the classic Lamb shift, and hyperfine 
splittings (hfs) in hydrogen, positronium, and muonium have 
been confirmed to better than 10 parts per million (ppm) (jll^h 
2 ppm yllj, 2 ppm C101> and 1 part in 10 million (j, respec- 
tively. The measurements of the electron and positron gyro- 
magnetic ratios agree with order-a 4 perturbative QED predic- 
tions to 1 part in 10 11 |4|. QED and gauge theory have thus 
been validated to extraordinary precision. 

In this paper we shall show how one can combine preci- 
sion atomic physics measurements to determine a fundamen- 
tal property of the proton to remarkable precision. The dif- 
ference between the ground-state hfs of hydrogen and muo- 
nium, after correcting for the different magnetic moments of 
the muon and the proton and for reduced mass effects, is due 
to the structure of the proton. The QED contributions for a 
pointlike nucleus essentially cancel. The largest proton struc- 
ture contribution to the hfs difference is proportional to the 
Zemach radius |5, 6], which can be computed as an integral 
over the product of the elastic electric and magnetic form fac- 
tors of the proton. The remaining proton structure corrections, 
the polarization contribution |7l |9l 1 1 Of! from inelastic 
states in the spin-dependent virtual Compton amplitude and 
the proton size dependence of the relativistic recoil correc- 
tions 

MM 

, have small uncertainties. As we shall show, the 
resulting high precision determination of the Zemach radius 
from the atomic physics measurements provides an important 
constraint on fits to accelerator measurements of the proton 
electric and magnetic form factors. 

An important motivation for examining form factor con- 
straints comes from the recent polarization transfer measure- 
ments of the proton electric form factor Ge(Q 2 ) dHHHIS 
The polarization transfer results are at variance with the pub- 
lished Rosenbluth measurements of Ge- The difference 
may well be due to corrections from hard two-photon ex- 
change fir)! [nil . One wants to examine with the maxi- 
mum possible precision whether the new determinations of 
Ge(Q 2 ), falling with respect to Gm{Q 2 ), is compatible with 



other information on the form factor. The extraction of the 
Zemach radius to be described here provides such a constraint. 

A sum rule for proton structure. We now show how one 
can use the hfs of the muonium atom (e~/i + ) to expose the 
hadronic structure contributions to the hydrogen hfs. For an 
electron bound to a positively charged particle of mass mjv, 
magnetic moment /ijv = (ffjv/2)(e/2mjv), and Lande g- 
factor gN, the leading term in the hfs is the Fermi energy, 
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Here, "N" stands for either the p or fi + nucleus. By conven- 
tion, the exact magnetic moment ^jv is used for the proton or 
muon, but only the lowest order term, the Bohr magneton ps, 
is inserted for the e~. 

The ground-state hydrogen hfs can be written as 



E h fs(e p) 



(1 + A QED + A p r + A s (2) 

+ ^hvp + ^pvp + ^weak)^f 



where Aqed represents QED corrections, A^. represents re- 
coil effects, including finite-size recoil corrections, As repre- 
sents the proton structure contributions, and A p hvp , A^ KJ) , and 
A^ cak represent the effects of hadronic vacuum polarization, 
muonic vacuum polarization, and weak interactions, respec- 
tively. The corresponding quantity for muonium is simply 



E his (e-[i + ) = (1 + Aqed + A 



^hvp 
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We define the fractional difference between the hydrogen 
and rescaled muonium hfs as 



»hfs 
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E Ms (e p) ^ (l + m e /m p ) 3 
EMs{e~p + ) p p (1 + m e /m M ) 3 

E his {e~ P )/E p F 
E Ms (e~p+)/E F ■ 



The large contributions from QED corrections cancel in Ahf s - 
Since the hfs of hydrogen and muonium, as well as the ratio of 
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muon and proton magnetic moments, have been measured to 
better than 30 ppb, A^fe can be determined to high precision 
from experiment. 

From Eqs. 10 and l|3}, we have 
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(1 + A QED + A P R + A s + A p hvp + A^ vp + A^ cak ) 
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Thus we can obtain a result for the proton structure contri- 
bution in terms of quantities measurable to high precision in 
atomic physics: 



Ac = A 



hfs 
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The cross terms are smaller than the uncertainties in the lead- 
ing terms, and here Aqed can be approximated as a/2n. 

The proton structure contributions consist of the classic 
Zemach term computed from a convolution of elastic form 
factors and the polarization contribution from the inelas- 
tic hadronic states contributing to the spin-dependent virtual 
Compton scattering: Ag = Az + A po \. In addition, as 
we discuss below, the relativistic recoil corrections of order 
am e /m p are modified by the finite size of the proton. The 
Zemach term takes into account the finite-size correction to 
the proton magnetic interactions as well as the finite-size dis- 
tortions of the electron's orbit in the hydrogen atom 0]: 
Az = — 2am e (r)z (l + S r z d ), where (r)z is the radius of 
the proton as calculated from the Zemach integral 
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with Ge and Gm the electric and magnetic form factors of 
the proton, normalized with G_e(0) = Gm(0)/(1 + k p ) = 
1, and k p = g p /2 — 1. Additionally, is a radia- 
tive correction to the Zemach term estimated in fiUl . It 
has been calculated analytically in |18] for the case where 
the form factors are represented by dipole forms: = 
(a/37r) [21n(A 2 /m 2 ) - 4111/420]. With A 2 = 0.71 GeV 2 , 
this yields <5™ d = 0.0153. 

The main part of the inelastic contribution can be con- 
structed from the work of Iddings [7] and Drell and Sulli- 
van 1 8]. Compact expressions are given by De Rafael |9], 
Gnadig and Kuti and Faustov and Marty nenko j^l in 
terms of the Pauli form factor F2 and spin-dependent struc- 
ture functions g\ and 52 of the proton. 

Evaluation of the constraint. We will consider each 
term on the right hand side of Eq. l|6). To compute 
Ahfs from we use the measured hydrogen hfs 11911 
E hl Je-p) = 1 420.405 751 766 7(9) MHz and muonium 
hfs I20I1 E his (e- fj,+) = 4 463.302 765(53) MHz. The mea- 
sured masses are m p = 938.272 029(80) MeV, m M = 



105.658 369(9) MeV, and m e = 0.510 998 918(44) MeV. 
The ratio of magnetic moments has been measured to high 
precision, ±0.028 ppm; the value obtained without input from 
the muonium hfs is 10 fip/fi p = 3.183 345 20(20). From 
these values we find Ahf s = 145.51(4) ppm. 

The recoil corrections A^ are separated into relativistic 
corrections A^ and additional radiative corrections A^ d . 
The order-a relativistic recoil correction A^ has been com- 



puted by Arnowitt 112 311 for muonium (N = /i). Bodwin and 
Yennie 1 1 1 ) quote the corrections to second order in a in their 
Eq. (1.10), which is analogous to Eq. l|8) below. Expres- 
sions for the radiative correction Af ad are given in J24I1 and 



12511 . With use of |21] a" 1 = 137.035 999 11(46) and J2€ 
= 0.001 165 920 8(6), the total correction is evaluated to 
be A£ = -178.34 ppm. 

Bodwin and Yennie full have also computed the correc- 
tions to their formula in the hydrogen case due to the finite 
size of the proton from elastic intermediate states. Note that 
these are finite-size corrections to the recoil correction and are 
distinct from the Zemach correction. A mark of the distinc- 
tion is that after scaling out the lowest order Fermi hfs, the 
recoil corrections go to zero as (m p /m e ) — > 00, whereas the 
Zemach correction does not. The Bodwin- Yennie pointlike 
result to order a 2 is El 
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with k p = 1.792 847 351(28). This gives Af cl = 
(—2.01 + 0.46) ppm, where the two terms are from O(a) 
and 0(a 2 ). Quoting 1 11], finite-size corrections change this 
to Af cl = (+5.22(1) + 0.46) ppm = 5.68(1) ppm, where 
the quoted error is an estimate using the dipole form fac- 
tor for the proton (both Ge and Gm) with mass parameter 
A 2 = 0.71 ± 0.02 GeV 2 . An additional radiative correc- 
tion 1 18] of 0.09 ppm brings A^. to 5.77 ppm. 

Volotka et al. 1 27] have reevaluated the finite-size correc- 
tions to the proton recoil corrections with the same magnetic 
radius, but with a charge radius taken from Ref. 1 28], and find 
A^, = 5.86 ppm, or 0. 18 ppm larger than Bodwin and Yennie. 
By forcing the magnetic form factor to reproduce their result 
for the Zemach integral, Volotka et al. obtain a second value 
of 6.01 ppm. We shall use the first Volotka result and include 
an uncertainty of 0.15 ppm to cover the difference between 
the modified Bodwin- Yennie and the second Volotka deter- 
minations. Note that structure-dependence uncertainty within 
the recoil corrections is still well under the uncertainty of the 
polarization terms, and that this uncertainty in the recoil term 
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TABLE I: Proton electric charge radius v/ Zemach contribution Az to the hfs, and Zemach radius (r)z for various parameterizations 
of Ge and Gm. The results should be compared to Az = —(40.0 ± 0.6) ppm or (r)z = (1.043 ± 0.016) fm, as obtained from analysis 

of atomic hfs data. The dipole form is Gm(Q 2 ) = (1 + k p )/(1 + Q 2 /0.71 GeV 2 ) 2 . The G E labeled JLab is 0] (i-0.13q^j) 
Parameterizations A-I and A-II are from Q. Those labeled Brodsky-Carlson-Hiller-Hwang (BCHH), I and II, use F 2 /F 1 = [1/k 2 + 
Q 2 /(1.25GeV) 2 ]" 1/2 and F 2 /F 1 = k p [1 + (Q 2 /0.791 GeV 2 ) 2 ln 71 (l + Q 2 /4m 2 )]/[l + (Q 2 /0.380 GeV 2 ) 3 ln 51 (l + Q 2 /4m 2 )], 
respectively j 30] . The last column gives the contribution to (r)z from Q > 0.8 GeV. 



Parameterizations 






A z 


(r)z (fm) 




G M 


Ge 


(fm) 


(ppm) 


total 


Q > 0.8 GeV 


dipole 


G m /(1 + k p ) 


0.811 


-39.32 


1.025 


0.310 


dipole 


JLab 


0.830 


-39.83 


1.038 


0.310 


A-I 


A-I 


0.868 


^tl.46 


1.081 


0.310 


A-I 


G M /{1 + Kp) 


0.863 


^41.34 


1.077 


0.310 


A-II 


A-II 


0.829 


^10.29 


1.050 


0.310 


A-II 


JLab 


0.855 


-40.70 


1.061 


0.310 


dipole 


BCHH-I 


0.789 


-38.85 


1.012 


0.310 


A-II 


BCHH-I 


0.816 


-39.72 


1.035 


0.309 


dipole 


BCHH-II 


0.881 


-40.91 


1.066 


0.310 


A-H 


BCHH-II 


0.905 


^11.77 


1.088 


0.310 



can be reduced as knowledge of the form factors improves. 

Estimates of the weak and vacuum polarization corrections 
are also given by Volotka et al. %Z~k . From these and from 
the individual values for Ahf s , A^,, and A P R , we obtain As = 
—38.62(16) ppm. Thus the contribution of proton structure is 
constrained by atomic physics with an uncertainty well under 
one percent. 

The Zemach term. We shall subtract the polarization con- 
tributions to isolate the Zemach term and then explore its 
relevance to new form factor parameterizations. Although 
the polarization contributions have been long known to be 
small |9j 1 1 (XI . the error in Az is essentially all due to the 
uncertainty in A po i. From Faustov and Martynenko fl, we 
take Ap i = 1.4 ppm ± 0.6 ppm, which implies Az = 
-(40.0 ± 0.6) ppm and thus (r) z = (1.043 ± 0.016) fm. 
The unit conversion used he = 197.326 968(17) MeV-fm. 

Predictions for Az and (r)z as computed from a selec- 
tion of parameterizations of the form factors are given in 
Table U The first row is the textbook standard, wherein 
both Gm and G e are given by the dipole form. The result, 
A z = -39.32 ppm = -38.72(1 + <5^ ad ) ppm, can already be 
found in full . New analytic fits to the form factors I29I l3dtl 
make a significant change in the Zemach integral, of up to 
6%. The form factor parameterization given in 1 28] yields |6] 
(r)z = 1.086(12) fm. It is not clear why the large difference 
exists. The scattering data is subject to radiative and other 
corrections; any difference highlights the usefulness of hav- 
ing the precise value that we have derived. Not all of the Az 
or (r)z for the different models in the table are compatible 
with the results extracted from the analysis of the atomic data. 
However, the Gm-Ge combination suggested in the third row 
from the end of the Table shows that fully compatible models 
exist. 



The table also shows results for the charge radius \J [r 2 E ) = 

'— 6^5-G 'e(Q 2 )\q 2 =o- The values compare to results from 
Lamb shift measurements [31] (0.871(12) fm), a continued- 
fraction fit to Ge S (0.895(18) fm), a standard empirical 
fit S (0.862(12) fm), and the 2002 Committee on Data for 
Science and Technology value (0.8750(68) fm). 

The differences among the Zemach integrals for different 
form factors derive mainly from the lower Q range of the in- 
tegral, where the different parameterizations of Ge are less 
variant, but their effect on the integral is greater. This is seen 
in the last two columns of Table[[] About 30% of the Zemach 
integral comes from Q above 0.8 GeV, but little of this has 
to do with the form factors. Recall that the numerator of the 
Zemach integrand is Ge Gm /(l + k p ) — 1, and for high Q the 
form factors fall away, leaving the "—1." In the region above 
0.8 GeV, the contributes 0.314 fm. 

Two fits by Arrington 1 29] are used in the Table, denoted 
A-I and A-II. Fit A-I uses only Rosenbluth data and A-II uses 
Ge/Gm from the polarization results 0. While 
A-II represents the data well overall, for Q below 0.8 GeV 
its Ge/Gm ratio falls too quickly by nearly a factor of two 
compared to the actual polarization data. The same is true for 
the fit denoted JLab 0. 

Discussion. In this paper we have shown how one can 
combine high-precision atomic physics measurements of the 
ground-state hydrogen and muonium hfs and the ratio of muon 
to proton magnetic moments to isolate the proton structure 
contributions. In our method, the theoretically complex QED 
corrections to the bound-state problem do not appear 13411 . 
Remarkably, the total proton structure contribution A5 = 
—38.62(16) ppm to the hydrogen hfs is determined to bet- 
ter than 1%. Since the polarization contribution can be de- 
termined from the measured spin-dependent proton structure 
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functions g\{x, Q 2 ) and g^ (x, Q 2 ), we obtain a precise value 
for the Zemach radius (r)z = (1.043 ± 0.016) fm, which is 
defined from a convolution of the Ge and Gm form factors. 
This new determination gives an important constraint on the 
analytic form and fits to the proton form factors at small Q 2 . 
The precision of the Zemach radius will be further improved 
when new, more precise data for gi and gi, especially at small 
v and Q 2 , becomes available. 

The proton structure terms can also be extracted using the 
hydrogen hfs alone 127113511 . The Zemach radius obtained this 
way is slightly smaller but consistent with our result. 

Conversely, by combining the muonium and hydrogen 
hfs data, one can obtain an "experimental" value for the 
purely QED bound-state radiative corrections: Aqed = 
1136.09(14) ppm. To minimize the uncertainty, we 
take advantage of the measured ratio l33ll m p /m e = 
1836.152 672 61(85). This value of Anm is consistent with 
the calculated QED correction used in l27U35ll . 

Our method of combining experimental atomic physics has 
other applications; for example, measurements of the differ- 
ence of the Lamb shifts (or Rydberg spectra) of muonium and 
hydrogen could potentially give a very precise value for the 
proton's electric charge radius, since again the QED radiative 
corrections essentially cancel. Similarly, the difference of lep- 
ton anomalous moments — a e directly exposes the hadronic 
and weak corrections to the muon moment. 
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